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Abstract 



1 Introduction 



> 

. We use an n-spin system with permutation symmetric zz-interaction for simulating 

I arbitrary pair-interaction Hamiltonians. The calculation of the required time overhead 

is mathematically equivalent to a separability problem of n-qubit density matrices. We 

derive lower and upper bounds in terms of chromatic index and the spectrum of the 

interaction graph. The complexity measure defined by such a computational model is 

related to gate complexity and a continuous complexity measure introduced in a former 

^ ■ paper. We use majorization of graph spectra for classifying Hamiltonians with respect 

^ , to their computational power. 

P . 

> 

^. 

, The most common models for quantum computers use single and two qubit gates as 

basic transformations in order to generate arbitrary unitary operations on the quantum 
registers. Most discussions about the generation of quantum algorithms, quantum 
codes and possible realizations had successfully been based on this concept. Mostly, 
even the definition of quantum complexity refers to such a model Nevertheless 
there is a priori no reason, why two qubit gates should be considered as basic operations 
for future quantum computers. In principle every quantum system could serve as a 
quantum register provided that its time evolution can be controlled in a universal way. 
At first sight, every definition of quantum complexity seems hence to be adequate only 
for a specific model of quantum computation. But it seems to be a rather general 
feature of Hamiltonians available in nature that particles interact with other particles 
in such a form, that the total Hamiltonian is a sum of pair-interactions. Therefore 
we want to base quantum complexity theory only on such a general feature.f] This 



^This might be seen in the spirit of D. Deutsch's statement "What computers can or cannot compute is 
determined by the laws of physics alone and not by pure mathematics." Chapter II 



1 



feature justifies the following control theoretic model: If n qubits are assumed to be 
physically represented by n particles, the only part of the system's Hamiltonian which 
can be changed by extern access is the free Hamiltonian of each qubit. These 1-particle 
Hamiltonians might be controllable since they are only effective Hamiltonians which 
are phenomenologically given by an interaction to many extern particles (mean-field 
approximation Q). Based on results of quantum control theory in 2-spin systems |Q 
we investigate the problem of simulating arbitrary pair-interaction Hamiltonians by a 
given one. We assume that the Hamiltonian of the n-system is a permutation invariant 
zz-interaction and show that the computational poweiQ of this Hamiltonian (together 
with local transformations on each spin) is at least as large as the power of quantum 
computers with 2 qubit gates. For infinitesimal time evolutions, it turns out to be even 
stronger. 

We develop a theory, where the computational power of a Hamiltonian for simu- 
lating arbitrary Hamiltonians is characterized by features of the interactions graphs. 
Standard concepts of graph theory like chromatic index and spectrum of the adjacency 
matrix together with majorization turn out to provide lower and upper bounds on the 
simulation overhead. Here we are interested in the exact overhead and not only in 
polynomial equivalence as in Q. 

2 Our model of computation 

Based on the approach of Q we consider the following model. The quantum system is 
a spin system, i.e. its Hilbert space is TCn = (C^)*^"", and its Hamiltonian G su(2"') 
consists only of pair- interactions, i.e. 

l<fc</<n 

where H^i acts only on the Hilbert space of the qubits k and /. We assume that for 
every k and / the Hamiltonian H^i describes a non-trivial coupling and is traceless. 
The system's Hamiltonian Hd is also called the drift Hamiltonian since it is always 
present. We assume that we can perform all unitaries in the control group K = 
SU{2) (8> • • • SU{2) arbitrarily fast compared to the time evolution of the internal 
couplings between the qubits. Let G be the unitary Lie group 5C/(2") and n G G be a 
unitary we want to realize. To achieve this all we can do is perform vi G K, wait ti, 
perform V2 G K, wait t2 ■ ■ , perform Vp (z K and wait tp. The resulting unitary is 

u = exp{iHdtp)vp ■ ■ ■ exp{iHdt2)v2 exp{iHdti)vi . 

This can be written as 

u = kp exp{ikj^Hdkp) ■ ■ ■ ex.p{iklHdk2) exp{ik\Hdki) 

where ki = Vi - ■ -vi for i = 1, . . . This is just the solution of a time-dependent 
Schrodinger equation with piecewise constant Hamiltonians - conjugates of the drift 

^in the sense of time required to generate unitaries 
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Hamiltonian H^i by unitaries of K - followed by the unitary kp E K. Let AdxiHci) 
denote the conjugacy class 

AdKiHd) = {Adk{Hd) = k^Hdk \k£K}. 



Definition 1 A continuous time algorithm A of running time T is a piecewise constant 
function t H[t) from the interval [0,T] onto the set AdxiHa) followed by some local 
unitary k ^ K . We say A implements u if u = ku[T) where (it(i))te[o,r] the solution 
of the time- dependent Schrddinger equation {d/dt)u{t) = —iH{t)u{t) with n(0) = /. 

The complexity of a unitary in this model is the running time of the optimal continuous 
time algorithms. 

Let denote the Pauli spin matrix cJq, (a = x, y, z) that acts on the ith spin. For 
simplicity, we assume that the drift Hamiltonian is 

l<fc<Z<?i 

The physical systems we have in mind might be for example solid states with long- 
range interactions. Of course one might object that the interaction strength always 
decreases with the distance between the interacting particles. It will turn out that 
the assumption on non-decreasing interaction strengths makes our model rather strong 
with respect to its computational power. One should understand our assumptions as 
the attempt to use a strong computational model which is still physically justificable. 
Many aspects of our theory can be developed in strong analogy for more general drift 
Hamiltonians. 

In Section |^ we will compare the computational power of our model with the power 
of quantum computers based on 2-qubit gates. Our arguments refer always to infinites- 
imal time evolutions, i.e., we will show that our model can implement quantum gates 
without overhead since we can simulate the time evolution implementing parallelized 
quantum gates. 

First we have to define what we mean by simulating the time evolution exp{iHt) 
during a small time interval [0, e] where H is an arbitrary pair-interaction Hamiltonian. 
Assume we have written H as a positive linear combination H = fJ-jHj with fij > 
and each Hj is an element of the conjugacy class AdK{Hd)- For small e the unitary 

JJexp(k^ji?j) 

i 

is a good approximation for 

exjp{ieH) = exp(ie^^ l^j^j) ■ 

3 

This approximation is implemented if the system evolves the time e/ij with respect 
to the Hamiltonian . The sum /i = [i^ is exactly the time overhead of the 
simulation. Hence the problem is to express as a positive linear combination such 
that the overhead [i is minimal. Of course such a procedure might not be optimal if one 
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were interested in the implementation of exp{iHs) for any special value of s. Here we 
want to imitate the whole dynamical time evolution (cxp(ii7s))s>o in arbitrary small 
steps e. Then the optimization reduces clearly to the convex problem stated above. 

In the following it will be convenient to use a concise representation for the drift 
Hamiltonian and the interaction to be simulated: a pair interaction Hamiltonian be- 
tween qubits k and I can be written as 



Hkl = ^ Jkl;a^O-a^p ■ 
a,l3=x,y,z 

The strengths of the components are represented by the pair-interaction matrix 



(3) 



Jkl 



Jkl;xx 'Jkl;xy 
Jkl;yx Jkl\yy 
Jkl;zx Jkl;zy 




p3x3 



(4) 



The total Hamiltonian H is represented by the J-matrix 



/ 



J 



J21 



J31 



\ Jnl 



Jl2 



•H2 



Jn2 



Jl3 



J23 



Jn3 



Jin \ 



J2n 



J3n 



])3nx3n 



(5) 



J 



To explain more explicitly, why our simulation problem is a convex optimization, 
we recall that every convex combination iJ,Hi -|- (1 — /x)iJ2 of two Hamiltonians Hi and 
H2 can be simulated with overhead 1 if Hi and H2 can. Remarkably, the problem of 
specifying the set of Hamiltonians which can be simulated with overhead 1 is related 
to the problem of generalizing Bell inequalities to ra-qubit states. More specifically, the 
convex problem can be reduced to the question 'how strong can 2-spin correlations be 
in a separable n-qubit quantum state?' 

Theorem 1 (Optimal simulation) The Hamiltonian H can be simulated with over- 
head fj, if and only if there is a separable quantum state p in (C^)®" such that 



A* 

where J denotes the J-matrix of H and I the 3n x 3n identity matrix. 

Proof: By rescaling the considered Hamiltonian, it is sufficient to show that this 
is true for all Hamiltonians in AdxiHa) with = 1. Assume we have written if as a 
convex combination H = fJ-jHj with Hj G AdxiHa). In order to show that there 
is a separable state of the desired form it is sufficient to show that the J-matrix of 
each Hamiltonian Hj satisfies the equation of the theorem for an appropriate separable 
state. Let Hj = uH^u^ . Hj can be represented by n three dimensional real unit 
vectors: to each qubit we associate the vector \Jk) = {Jk-x, Jk;y, Jk;zY S where 
Uk(^zu\ = Jk;x(^x + Jk;y'^y + Jk;z'^z and « = tii Cg) . . . (8>u„. The pair-interaction matrices 
are given by the matrix products Jki = \Jk){Jl\- 
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By the Bloch sphere representation we have a correspondence between the unit 
vectors \Jk) and the projections pk in defined by J^-^a = ^^{Pk'^a)- Let p be the 
product state p := pi . . . (g) pn- Then we have Jki-^ap = tr(p(T^(T^) for all k ^ I. Note 
that the product of two different Pauli matrices is the third Pauli matrix multiplied 
by a scalar. The only problem that remains is that we may have ti{pa^a^) ^ for 
a ^ f3. We substitute p by a state p in such a way that the expectation values of 
all traceless 1-qubit observables vanish and the expectation values of all considered 
2-qubit observables remain unchanged. For every |Jfc) we can find G SO (3) such 
that U'j^\Jk) = —\Jk)- This rotation corresponds to conjugation of the qubit A; by a 
unitary n'^. To —\Jk) corresponds the projection p'/, := I2 — Pk- Let 

p:=^{pi(^...(^pn + p'i(S>...(^p'J. 

Then we have tr(/)cr^(T^) = tr(/>(T^(T^) for all A; 7^ / (all vectors are multiplied by — 1 
and therefore there is no effect on the pairs) and Jkk is the 3x3 identity matrix. 

Assume conversely we have a separable state of the desired form. Take its decom- 
position into pure product states. By the Bloch sphere representation we obtain the 
required conjugations of the drift Hamiltonian Ha- The time they have to be applied 
are given by the coefficients in the convex decomposition. □ 

3 Lower and upper bounds 

A simple lower bound on the simulation time overhead can be derived from the fact 
that J/p + 1 has to be a positive matrix, which is an easy conclusion from Theorem |l[ 

Corollary 1 (Lower bound) The absolute value of the smallest eigenvalue of the 
J -matrix is a lower hound on the simulation overhead of H . 

Proof: The matrix 

(tr(pf7^a^))H;a/3 

is positive for every state p in (C^)®": let \d) = {dk;a) be an arbitrary vector and 
^ = Y,k,a dk;a(^a- Then we have 

dk;at^{pa'^a^f^)di,f3 = tT{pAA*) > 0. 

k,l,a,P 

□ 

Now we show that our computational model is at least as powerful as the usual 
model with 2-qubit gates, even if also cares about constant overhead. We describe here 
briefly the quantum circuit model and introduce the weighted depth following [p!l|]. It is 
a complexity measure for unitary transformations based on the quantum circuit model. 
We assume that two qubit gates acting on disjoint pairs of qubits can be implemented 
simultaneously and define: 

Definition 2 A quantum circuit A of depth k is a sequence of s steps {Ai, . . . ,As} 
where every step consists of a set of two qubit gates {uki}k,i acting on disjoint pairs 
(k, I) of qubits. Every step i defines a unitary operator Vi by taking the product of all 
corresponding unitaries in any order. The product u := ILi^gi'i is the 'unitary operator 
implemented by A '. 
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The following quantity measures the deviation of a unitary operator from the identity: 

Definition 3 The angle of an arbitrary unitary operator u € 5[/(4) is the smallest 
possible norrr^ \\a\\ of a self-adjoint operator a G 5u(4) which satisfies exp(m) = u. 

It coincides with the time required for the implementation of u if the norm of the used 
Hamiltonian is 1. We consider only the angle of two-qubit gates, i.e. we do not include 
the angle of local gates in the definition of the weighted depth. The notion of angle 
allows us to formulate a modification of the term 'depth' which will later turn out to 
be decisive in connecting complexity measures of discrete and continuous algorithms: 

Definition 4 Let Oi be the maximal angle of the unitaries performed in step i. Then 
the weighted depth is defined to be the sum a = ^ ^ Ui . 

Assuming that the implementation time of a unitary is proportional to its angle, 
the weighted depth is the running time of the algorithm. We first need two technical 
lemmas to show that such an algorithm can be simulated by our computational model 
with complete zz-Hamiltonian without any time overhead. 

Lemma 1 Let M be a set of qubit pairs, such that no two pairs contain a common 
qubit. Then we can simulate 

hm= 

with overhead 1. 

Proof: This has been noted in [Q. Theorem ^ proves a more general statement. □ 

Lemma 2 Let = (Tz be the drift Hamiltonian of a 2-spin system. All Hamil- 
tonians H € 5u(4) can be simulated with overhead less than \\H\\. 

Proof: We first assume that H contains no local terms, i.e. H = ^ JapCTa ® C/j. 
Let Ji2 be the matrix representing H. Conjugation oi H hy k = w^v £ SU(2)0SU{2) 
corresponds to multiplication of Ju hy U € SO (3) from the left and hyV G 50(3) from 
the right. By the singular value decomposition there are U,V £ SO{3) such that 
Ji2 = Udiag{sx, Sy, Sz)V where Sx,Sy,Sz are the singular values of J12. Equivalently, 
there is A; G SU{2) (g) SU{2) such that kHk"^ = Hs^,sy,s^ where Hs^^sy,s^ = s^Cx ^ (Tx + 
Syay®ay + SzOz®(Tz- By computing the eigenvalues we see that ||-ffs2;,sy,sz II = Yla \^a\- 
The simulation time overhead can not be more than the right hand side since each 
term SaOa ® era can be simulated with overhead Sq,. 

Let H contain local terms, i.e. H = JaaC^a ^Ca + lf^a + fe^l. We can split 
H = H' + H" where H' is the non-local part and H" the local one. By the Trotter 
formula we can simulate the parts independently. The simulation of H" takes no time 
by assumption. It remains to show that ||-ff'|| < ||-ff||. We may assume that H is 
invariant with respect to qubit permutation since ||^-ff + ^-ffexll ^ ||^|| where H^x is 
the Hamiltonian obtained from H by exchanging the qubits. By conjugation we can 
obtain a Hamiltonian of the form H = Hs^^sy,s^ + ^(l ax + cTx ^ 1). By computing 
the eigenvalues we see that ||-ffs^,sy,s^ || < ||-ff||- D 

^Here ||.|| denotes the operator norm given by ||a|| := maxx \\ax\\ where x runs over the unit vectors of 
the corresponding Hilbert space. 
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Corollary 2 Let A = {uki}k,i be a step of a quantum circuit and a its weighted depth. 
Then the zz-model can simulate the unitary implemented by A with overhead a. 



Proof: Let M = {{k,l)} be the set of the pairs which the two-qubit gates act 
on. No two pairs in M contain a common vertex and therefore we can simulate the 
Hamiltonian Hm = '^(k i)gm "^z^z with overhead 1. Let H^i be the Hamiltonian of 
minimal norm such that u^i = exp{iHki) for every {k,l) G M. Now we can simulate 
every H^i parallely with overhead less than \\Hki\\ by conjugating Hm- ^ 

Our goal is to compare interactions with respect to the simulation complexity in 
our model given by the complete zz-interaction and the quantum circuit model. For 



doing so, we need some basic concepts of graph theory [|10[. A graph is an ordered 
pair G = {V, E) with V C {1, 2, . . . , n} and = {ei, 62, . . . , Cm} 'xV. Elements 
of V are called vertices. They label the qubits. Elements of E are called edges. They 
label the pair-interactions between the qubits. An edge e = {k, I) is an ordered pair 
of vertices k and I called the ends of e. We consider only undirected graphs with no 
loops. To have a unique representation we require that k < I. Two distinct edges are 
called adjacent if and only if they have a common end vertex. A subset M of the edge 
set E is called independent if no two edges of M are adjacent in G. A graph G is 
called complete if every pair of distinct vertices of G are adjacent in G; such a graph 
is denoted by Kn- Rephrased in this language, our drift Hamiltonian is of the form 

{k,l)GE{K„) 

and is called in the following the complete zz-Hamiltonian. 

Definition 5 Let H be an arbitrary pair-interaction Hamiltonian. For every non- 
negative real number r we define the interaction graph Gr as follows: Let the qubits 
label the vertices and let the edges be all the pairs {k,l) with the property 

\\Hk,l\\ > r. 

The chromatic index x' is the minimum number of colors permitting an edge-coloring 
such that no two adjacent edges receive the same color or equivalently a partition 
E = Ml U M2 U . . . U M^i into independent subsets of E. The following quantity turns 
out to be an upper bound on the overhead. 

Definition 6 We define the weighted chromatic index of H 

POO 

X' ■■= / Xrdr (7) 







where x'r denotes the chromatic index of Gr- 
in a former paper we have introduced the weighted chromatic index as a com- 
plexity measure of the interaction. This point of view has been justified by two argu- 
ments, where the first one is an observation in [11 1: 



Theorem 2 The evolution generated by a pair-interaction Hamiltonian H during the 
infinitesimal time period dt can be simulated by a parallelized 2-qubit gate network with 
weighted depth x' dt if x' the weighted chromatic index of H . 
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The second argument to consider chromatic index as a complexity measure for the 
interaction is only intuitive: in general, it should be easy to control interactions on 
disjoint qubit pairs, whereas one should expect that its unlikely that one can control 
simultaneously the interaction between qubit 1 and 2 and the interaction 1 and 3 at 
the same moment. This 'a priori'-assumption of |jlT| can be partly justified by the 
following corollary which is an easy conclusion of Corollary |2| and Theorem ^. 

Corollary 3 The time overhead for simulating the Hamiltonian H in the zz-model is 
at most the weighted chromatic index of H . 

The assumption that the drift Hamiltonian contains only pair-interaction of the 
form (Tz^CTz can be dropped. Let H = J2a p Jcip^a®(^p be an arbitrary pair-interaction. 
By conjugating H with {I®/, /(gjcx^, o"^ (8)/, (^ctz} we obtain JzzO'z^f^z- This can be 
done with overhead 1. The bounds of the corollary must be divided by the minimum 
Jzz of all pair-interactions occurring in H. 

4 Applications 

The graph theoretical nature of our optimization problems becomes even stronger if 
we reduce our attention to one type of interactions, namely zz-interactions. Then the 
desired Hamiltonian is completely described by a weighted graph. 
We consider the problem to simulate the time evolution 

H = Y,JkicTWz- (8) 
ik,l) 

when the complete zz-Hamiltonian is present. We first show that in this case it is 
sufficient to use conjugation by ax only. Let H' := az ^cTz- Note that {ax ®I)H'{ax (8) 
/) = —H' and {ax ® ax)H'{ax <8> ax) = H' . In the following we denote conjugation 
by ax by — and no conjugation by -|-. The Hamiltonian to be simulated contains only 
terms of the form Jki;zz(^z'^i by assumption. If it is written as a convex combination 
of elements of AdK{Hd) it is sufficient to show that for each of these elements there 
is a procedure which cancels the terms Jki-af3 for (a, /?) / {z, z) without any effect 
on the Jki,zz terms. Therefore consider H = kH^kK Then we can also achieve the 
Hamiltonian Hzz = Yl{k i) Jki;zz(^z'^i with overhead 1. For every qubit i there is a Ji-z 
such that Jki-zz = Jk;zJi;z for all edges {k,l). We express each Ji-z = cf — c~ with 
< c^, < 1 and + = 1- Let K = {I, ax} ^ . . . ® {I, ax}- We conjugate the 
drift Hamiltonian hy u = ui ® U2 ® ■ ■ ■ ® Un € K for time t{u) = ]Xi=i Ci{u) where 
Ci{u) = c'^ \i Ui = I and Ci{u) = if Uk = ax- We have 

^ t{u)uHu^ = Hzz ■ 

u&K 

Since we restrict our attention to interactions with zz-terms only a shorter notation 
will be useful. To each edge e = (/c, /) of G, we associate a real number wu called the 
weight of e. The resulting graph is called a weighted graph. Its adjacency matrix J is 
the real symmetric matrix with zeros on the diagonal defined by 

Jii := , Jfci := Wki and Jik ■= Wki (9) 



8 



for all edges {k, I) of G. An unweighted graph can be considered as a weighted whose 
edges all have the weight 1. 

A (unweighted) graph is bipartite if its vertex set can be partitioned into two 
nonempty subsets X and Y such that each edge of G has one end in X and the other 
in Y. The pair (X, Y) is called a bipartition of the bipartite graph. The complete 
bipartite graph with bipartition {X,Y) is denoted by G{X,Y). 

A Seidel matrix defines a modified adjacency matrix S = {ski) for (unweighted) 



graphs in the following way [12]: 



Ski 



-1 if A; and / are adjacent k ^ I 
1 if A; and / are non- adjacent 



and Sfcfc = 0. Obviously, S = K — I — 2 J, where K denotes a square matrix all of whose 
entries are equal to 1 and J the adjacency matrix of G. 

Theorem 3 (Optimal simulation) A graph G can he simulated with overhead 1 if 
and only if it can he expressed as a convex combination 



J = ^tA (10) 



where the sum runs over the Seidel adjacency matrices of all complete bipartite graphs, 
i.e., over 2"'~^ possible matrices. 

Proof: By assigning to each vertex either + or — we have a bipartition of the vertex 
set: X contains all vertices with + and Y all vertices with — . The sign of the edge 
{k, I) is — if and only if the edge has one end in X and the other end in Y and + 
otherwise. The edges with — define the complete bipartite graph G{X,Y). We also 
include the case X = and Y = V to cover the case when + is assigned to all knots. 
Therefore all we can achieve in the one step is K — I — 2J(X, Y). □ 

Corollary 4 (Lower bound) The absolute value of the smallest eigenvalue of the 
J -matrix is a lower bound on the simulation overhead. 

We present now some upper bounds on the overhead. A graph G = (y',E') is 
called a subgraph of G' if y C y and E' C E. A clique of G is a complete subgraph 
of G. A clique of G is called a maximal clique of G if it is not properly contained in 
another clique of G. A clique partition P of G is a partition of E[G) such that its 
classes induce maximal cliques of G. Given a set C oi h colors, an /i-coloring of P in 
G is a mapping from P to C, such that cliques sharing a vertex have different colors. 
Let the clique coloring index c(G) be the smallest h such that there is a partition P 
permitting an /i-coloring ||l^. We say the graph G consists of independent cliques if 
c(G) = 1. 

Lemma 3 (Upper bound) Let G he a graph consisting of independent cliques. We 
can simulate the Hamiltonian He with overhead 1 which is optimal. 
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Proof: Let w > 2 be the number of maximal cliques. We construct uj vectors Sj of 
length 2'^~^ as follows: 

si = (+ + + + + + + + •••) 

52 = (+- + - + - + ) 

53 = (+ + -- + + ) 



(+ + 




where + stands for 1 and — for —1. The scalar products are {si,Sj) = 2^~^6ij. We 
partition the time interval into 2'^^^ intervals of equal length. In the mth interval we 
conjugate all qubits of the ith clique by ax if Sj^m = — and do nothing otherwise. This 
is optimal since q < —1 where q is the smallest eigenvalue of G. □ 
The scheme used in the proof is time optimal. But the number of conjugations 
grows exponentially with the number of cliques w. It is possible to use the conjugations 
schemes based on Hadamard matrices [^. There the number of conjugations grows only 
quadratically with uj. 

Corollary 5 (Upper bound) Let G be an arbitrary graph. Then the Hamiltonian 
Hg can be simulated with the overhead c[G). 

Note that if M is an independent set than the graph G = {V, M) consists of independent 
cliques. Therefore the chromatic index is an upper bound on clique index. However, 
this bounds is not always good. Consider e.g. the graph G containing all edges that 
do not have 1 as end vertex. Then the chromatic index of G is still high but the clique 
coloring index is only 1. 

Since the optimal simulation of graph consisting of independent cliques has overhead 
1 one might think that the clique index is the smallest overhead. But this is not so as 
shows the following example. Consider the star G = {V, E) with V = {1, . . . , 5} and 
E = {(1, 2), (1, 3), (1, 4), (1, 5)}. The clique index of is 4 but the optimal simulation has 
overhead 2 only. The vectors can be chosen as si = (+ + ++), S2 = ( — h ++), S3 = 

(H S4 = (+ H h), S5 = (+ + H — ) and each of the four intervals has length 

1/2. This is optimal since the smallest eigenvalue of the adjacency matrix of G is —2. 



5 Quasi-order of Hamiltonians 

Let H and H be arbitrary pair-interaction Hamiltonians. We investigate the question 
whether H can be simulated by H with overhead Note that this defines a quasi- 
order of the pair- interaction Hamiltonians for = 1. A partial characterization of the 
quasi-order is expressed in terms of majorization of the spectra of the corresponding 
matrices J and J. Similar methods have been used to derive conditions for a class of 
entanglement transformations and to characterize mixing and measurement in quantum 
mechanics |^ |^ . 

Suppose that x = (xi, . . . , Xd) and y = (yi, . . . , yd) are two dimensional real vectors. 
We introduce the notation [ to denote the components of a vector rearranged into non- 
increasing order, so x'^ = (x|, . . . where (x| > X2 > . . . > x^^). We say that x is 
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majorized by y and write x ~< y, if 



i=i i=i 

for k = 1, . . . ,d — 1, and with equahty when A: = d Q. 

Let Spec(X) denote the spectrum of the hermitian matrix X, i.e. the vector of 
eigenvalues, and X{X) denote the vector of components of Spec(X) arranged so they 
appear in non-increasing order. Ky Fan's maximum principle Q states that for any 
Hermitian matrix A, the sum of the k largest eigenvalues of A is the maximum value 
tr (AP), where the maximum is taken over all fc-dimensional projections P, 

k 

^ Xj(A) = maxtr(AP) . 
i=i ^ 

It gives rise to a useful constraint on the eigenvalues of a sum of two Hermitian matrices 
C := A + B, that A(C) -< X{A) + X(B). Choose a fc-dimensional projection P such that 

k k k 

^ X,{C) = tr(CP) = tr(^P) + tT{BP) < ^ A,(^) + X,{B) . (11) 
j=i j=i j=i 

This permits us to derive a lower bound on the simulation overhead. 

Lemma 4 (Majorization) Let H and H he arbitrary pair-interaction Hamiltonians. 
A necessary condition that H can be simulated with overhead fj, by H is that Spec( J) ~< 
/iSpec( J). 

Proof: By representing the Hamiltonians by their J-matrices we see that H can 
be simulated with overhead // if and only if there is a sequence of orthogonal matrices 
Uj = Uji® ... Ujn G SO{3) © ... © SO{3) and fij > with ^j = such that 



The proof now follows from the inequality (11). 



We consider now the problem to reverse the time evolution ex.p{iHdt), i.e. what is 
the overhead of simulating —H^i when is present. 

Lemma 5 (Lower bound on inverting) Let r be the greatest eigenvalue and q the 
smallest eigenvalue of J. Then fj, > ^ is a lower bound on the overhead for simulating 
-Hd by Hd. 

Proof: This is a direct consequence of the Weyl inequality (see Theorem III.2) 
Xd{A + B) > Xd{A) + X[B) for the sum of two Hermitian matrices where X^ denotes 
the smallest eigenvalue. □ 
Let G be a connected graph and Hd = '}2{ki)&E{G)'^^'^^z- ^ is not connected 
then the components can be treated independently. For the spectrum the following 
statements hold (see ||l^, Theorem 0.13): 

1 < r < n — 1 , — r < q < —1 . 



11 



It is interesting to note that this gives a tight bound for simulating —H^ when G = Kn 
since for a complete graph we have r = (n — 1) and q = —1. An upper bound is the 
(weighted) chromatic index x'iKn) which is either n (if n is even) or n — 1 (otherwise). 
This simple example shows that the inverse of the natural time evolution may have a 
relatively high complexity. 

Lemma 6 Let the drift Hamiltonian Hd he an arbitrary pair-interaction Hamiltonian. 
If the interaction graph Gq{H(i) is bipartite then we can invert the time evolution with 
overhead of less than 3. 

Proof: Let S = {cTrc,ay,az}- We have Ylu^s'^^''^^ ~ ~^ ^"-"^ ^ ^ su(2). Let X,Y 
be the bipartition of Go{Hfi). By conjugating all qubits in X with elements of S we 
obtain —Hd. □ 

If H contains only cr^ ^ Oz then the overhead is 1. This is optimal since we have 
^ > 1. 
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